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LIBRARY ROUTINE S 3 - 130 



Logarithm (DOI or SADOl) 

Closed subroutine with standard entry 

14 

°> 1, 2 

-38 
Jfeximum error 2 ^ assuming exact argument * Has 

-39 
biased error from + to -2 o 

Approximately 51 milliseconds 

This routine is entered with the argument x in A, 

x being within the limits < x < 1. Log p x is 

left in AQ with the integer part in A and the 

fractional part in Q with q^ = I » If x = the routine 

will loop for determining the integer part of the logarithm 

and thus if it is possible that x - this should be 

tested for before entering the routine • 

This result may be converted to any desired base 

with the use of the formula log_x = log^x log 2. The 

"a ~c a 

two most often used conversion factors ares 



log 2 = 0.69314 71805 60 
log 1Q 2 = 0.30102 99956 64 

Since log p x is left in AQ, a right shift may be necessary 
to bring the integer part of the logarithm into Q before 
multiplying by the conversion factor- In this process 
q_. remains equal to one and thus the logarithm has the 
correct sign. In general* for < x < 1* it will be 
sufficient to shift right 6 times since the integer part ' 
can be no less than -39* 

The direct conversion factors after shifting right 
six places for base e and base 10 respectively ares 

(64/100) log 2 = 0,44561 41955 58 

(64/100) log 1Q 2 = .,19265 91972 25 

thus leaving (l/lOO) log x or (l/lOO) log in x in the 



'10* 



accumulator after multiplying . 



- 2 - 

The two most significant decimal digits are interpreted 
as the integral part of the logarithm- 
If l/2 < x < 1, the correct result is in Q alone when 
leaving the routine and no shifting is necessary before 
base conversion. More accuracy can thus he retained 
in the final result „ 
MATHEMATICAL METHOD The quantity log 2 x is found in this routine to take 

advantage of the binary operation of the computer . 
The integer part of log p x is found by considering 
x as 2 (w) where < x < 1, 1/2 < w < 1, and m is 
a. positive integer . Then 

log 2 x = 10g 2 2~ m (w) = -m + log ¥ 

log w is fractional and -m is the integer part of logpX. 

The quantity w is then used to calculate the fractional 

part of log p x from the series 

39 -i 
log^w = -1 .+ .51 a. 2 where a. is either 

zero or one * 

The quantities a, are determined from the recurrence 

relation 

p = x 



p i + i ■ *± * *i * x ' z 

= 2P? if F? < 1/2, 



1 1 

The a. are found one at a time by the relations 

a... = 1 if if > 1/2 



*i+l ~ " A x i 



a i+1 = if 



F? < 1/2 
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0QK(S3) 




1 





40 F 


Store x in location 






K5 F 


- 




1 


42 13L 
41 2F 


Plant link 




2 


50 F 








22 5L 






3 


40 F 








L5 2F 


Form "0" digit for fractional part 




4 


26 11L 








00 IF 






5 


F5 IF 


- 






40 IF 


J Count for integer part 




6 


SJ F 


IT 






32 4L 


Test for > l/2 




7 


S5 F 
2b 9L 






8 


50 F 


_ 






7J F 


J Form r 




9 


40 F 








L4 F 






10 


36 3L 


Test for P7 > l/2 






F5 2F 


Form n l M digit for fractional part 




ll 


L4 2F 
40 2F 






12 


36 8L 








50 2F 


Put fractional part in Q 




13 


Fl IF 


Put integer part in A 






22 F 


Exit 





